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Abstract 

Random jitter in the soliton arrival time (the Gordon-Haus effect) is analyzed for solitons being solutions 
of the integrable modified nonlinear Schrodinger equation. It is shown that the mean square fluctuation 
of the soliton position depends on the soliton parameters which can be properly adjusted to suppress the 
Gordon-Haus jitter. 

1 Introduction 

Long-distance soliton-based fibre transmission meets with difficulties imposed by random jitter in the soliton 
arrival time caused by the spontaneous amplifier noise (the so-called Gordon-Haus effect |l|]).The amplifier noise 
incorporated by the soliton produces a random soliton frequency shift which leads to a timing shift because 
of the group velocity dispersion. The Gordon-Haus effect limits error-free propagation of the soliton. There 
is a few ways to partially overcome the Gordon-Haus limit, including the use of linear filtering |J, the 
dispersion compensation means ^, |6) , the use of a sequence of two different media to reduce a path-averaged 
dispersion (7[ [|. 

From the viewpoint of the inverse scattering transform, the theory of the Gordon-Haus effect is a direct 
consequence of the adiabatic soliton perturbation theory for the nonlinear Schrodinger (NLS) equation, with 
perturbation being the amplifier noise ||. The NLS equation serves as the integrable model describing the 
picosecond soliton dynamics in fibres. On the other hand, when dealing with ultrashort optical pulses with 
duration < 100 fs, the NLS equation should be modified to adopt more subtle effects, such as the nonlinearity 
dispersion, the Raman self-frequency shift and the third-order dispersion |Io[| . It is remarkable that the account 
for the nonlinearity dispersion does not break the integrability of the equation. In other words, the modified 
NLS (MNLS) equation 

iu z + ^u tt + \u\ 2 u + ia{\u\ 2 u) t = 0. (1) 

is still integrable by means of the inverse scattering transform p"T[ , though the linear spectral problem associated 
with the MNLS equation differs from that for the NLS equation (with a = 0). Here u is the normalized slowly 
varying amplitude of the electric field envelope, z and t are the normalized propagation distance and time 
in the frame comoving with the group velocity, the real parameter a governs the effect of the nonlinearity 
dispersion. We consider the MNLS equation (|lj) as the integrable model for ultrashort optical pulses, i.e., 
playing the same role as the NLS equation does for picosecond solitons. Thereby, we change the status of the 
nonlinearity dispersion term from being a perturbation in the NLS equation to the essential ingredient of the 
MNLS equation. It is important that such a change is in no way an issue of our convenience. It was shown by 
Ohkuma et al. fl2]| that numerical simulation of soliton propagation revealed a number of features which cannot 
be accounted for by treating this term as a perturbation term in the NLS equation. In the recent paper [jl3| , the 
adiabatic perturbation theory for MNLS solitons was elaborated. So, a question to analyze the Gordon-Haus 
effect for the MNLS solitons arises naturally. 
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In this communication, we derive analytically the mean-square displacement fluctuation for the MNLS 
soliton propagating in a fibre. It follows from our results that the Gordon-Haus effect for the MNLS solitons, as 
distinct from the NLS solitons, can be suppressed by properly adjusting the soliton/fibre parameters, without 
making use of external means. In the end of the paper we briefly discuss the suitability of this result to actual 
femtosecond solitons. 



2 Formulation of the problem 

We consider the perturbed MNLS equation 

iu z + -Utt + \u\ 2 u + ia(\u\ 2 u) t = s(z, i), (2) 

where s{z, t) stands for spontaneous amplifier noise. In order to digress the details of minor importance, we 
consider distributed gain that exactly compensates for the fibre loss. Further, we consider the amplifier band- 
width to be much larger than the soliton bandwidth. Because erbium-doped amplifiers have a gain bandwidth 
of the order of 40 nm |5), this assumption works well for the ultrashort pulses. We consider the noise s(z,t) as 
the delta-correlated function both in time and space, 

< s(z, t)s(z', t') >= A5(t - t')8{z - z'). (3) 

The estimation of the coefficient A can be found in Q. Finally, we treat the noise s{z,t) as being small to 
justify the perturbative approach. 

The unperturbed soliton of the MNLS equation ([!]) was derived as early as in 1983 by Gerdjikov and 
Ivanov jt4j. We will use more simple and transparent expression 1 1 3 for the soliton solution of (|l|) : 

w (ke x + ke x ) 1 

Here x and ip are linearly expressed through coordinates z and t : 

x = -—+q(z), ip = vt + <f)(z), q(z) = a+—z, <f>{z) = Lp - \{v 2 - —z)z. (5) 

Parameters a and ip determine initial position and initial phase of the soliton, the complex parameter k = £ — ir), 
£, 1] > determines velocity (more exactly, shift of the reciprocal soliton velocity) v and width w of the soliton, 

v(z) = ±-*(k 2 + P), w(z)=- J a (6) 
a a 2(k z ~ k 2 -) 

where the z-dependence of the soliton velocity and width arises from the perturbation-induced z-dependence of 
the parameter k (see below (0)). Functions q(z) and <p(z) are the direct analogs of the Gordon's position and 
phase parameters for the NLS soliton p5| . 

The soliton ([|) has a number of peculiarities which distinguish it from the NLS soliton jt3|]. For example, 
the important invariant of eq. (pj), namely, the optical energy 
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has the upper limit 27r/a. It should be stressed that a enters the denominator of the soliton solution (|4|) (because 
of w) , so we account non-perturbatively for the nonlinearity dispersion. In other words, the parameter a does 
not to be small in general. Nevertheless, a non-trivial limiting procedure exists p3| permitting to restore 
from (Q) the NLS soliton as a — > 0. Indeed, let us take the limit 



k = \- f fc NLS + 0{a 2 ) 



&NLS = £NLS + iVNLS- 



(7) 



Then it is easy to see that the soliton velocity and width (j^) transform in the limit (Q) as v — » 2£nls, w — > 
1/2t7nls, i- e., to the NLS parameters, while the MNLS soliton ([|) produces exactly the NLS soliton unls = 
2i7?NLsexp(i'0NLs)sechx N Ls with x NLS = a - 2t]- n - ls (t - 2£ NLS 2:), ^nls = <t> + 2£nlst - 2(^ LS - ??nls) z - 

As it was shown in ||l3| , a perturbation- induced z-evolution of the parameter k is given by the simple formula: 



dk i 2 
— — = —ak 

dz 2 



dx- 



s(x)e-^ + «(-z)e* v ] 



(fee 1 + ke x ) 2 
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(8) 



which is transformed in accordance with (o) into the corresponding relations for velocity and width: 
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dx 
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Here S v (z) and S w (z) stand for the noise sources driving velocity and width. Because noise s(z,t) is given in 
terms of the correlation function (^|), the responses are also expressed in the same form. 

We are mostly interested in the mean-square fluctuation (q(L)q(L)) (L being a fibre length) of the soliton 
displacement q(z) (j^) which is given in the presence of perturbation by 



q(z)=a + / dz'v(z )w~ L (z r ) 



with v and w determined from (|9|). The jitter of the soliton displacement in the comoving frame is calculated 
from 

? = £ + (10) 

dz dz w 

where z-evolution of a is given by fll3|] 



da da+ da- _ da + 

T = —i ^-j— = ^a(z), — — = wq(z)b w {z), 

dz dz dz dz 



(11) 
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Since the soliton velocity and width enter eq. ( O ) in the combination vw 1 , we find from (|^) the noise source 

oo 

T-(-) = -2k 3 (v + -) [ dx- ^ - X9 [s{x)e-^ + s(-x)e^} 

dz y w J w' J (ke- x + ke x ) 2 

— OO 
OO 

_ 2 F(w--) / dx - _l* [s(x)e-^ + s(-x)e^] = S vw -i(z). (12) 
w J (ke x + ke x y 

— oo 

Finally, note that the mean-square fluctuation of the quantity B caused by the noise source Sb is determined 

by 

z z' 

(B(z)B(z')) = ( J dCS B (C) J dCS B (C)). 



Now we have everything to calculate the mean-square fluctuation of the soliton displacement. 

3 Results 

It is seen from (|l^) that the mean-square fluctuation of the soliton displacement is produced by both the noise 
source S a driving the displacement directly and the fluctuations of the combination vw~ 1 . We consider noise 
sources as independent, so the fluctuations they produce are additive. This gives 

(q(L)q(L)) = (q(L)q{L)) a + (q(L)q(L)) vw -i . 

Here 

L L 

{q(L)q(L)) vw -i = fdz f dz'(-(z)-(z')). (13) 

J J WW 

o o 

and 

L L 

(q(L)q(L)) a = J dz j dz'(S a (z)S a (z')). 



In turn, the correlation function for the velocity-to-width ratio is determined by the noise source S vw -i (12): 

z z' 

{-{z)-(z')) = dC [ d('(S vw -,(OS vw -i(t')). (14) 

WW J J 



Straightforward calculation by means of eqs. (§) and © gives 

(S vw -i (z)S vw -i (z')> = AF(a, v, j)5{z - z') (15) 

with 

4 

F(a,v,j) — — (1 — av){(l — 2oev )[(1 - aw) sec 2 27 - 2(1 - lav)\ 

+ ([(1 - aw) sec 2 2 7 - 2(1 - 2av)} 2 + (1 - 2ai>) 2 sec 2 27)7 cot 27}. (16) 
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Figure 1: Profile of the function F(a, V, 7) determining the correlation function for the soliton velocity-to-width 
ratio for a = 0.05. The soliton width w is related with the parameter 7 through (|l7|). 



In the course of calculation, we ignore, because of smallness of noise sources, small variations of the parameter 
k. Thereby, the parameters V, w and 7 refer to their initial values. (Remind that distributed gain compensates 
exactly for fiber losses). Besides, we express the soliton width w in terms of v and 7 due to the relation 



Finally, we obtain from ( |14| ) and (|T^) the correlation function for the soliton displacement driven by fluctuations 
of the velocity-to- width ratio: 



Similar calculation for the ^-driven correlation (q(L)q(L)) a in accordance with ( ]TT| ) displays the linear growth 
with L. Hence, for large propagation distances the correlation function ([l8|) provides the main contribution to 
the total mean-square fluctuation of the soliton displacement. 

4 Discussion 

The above dependences on L of the mean-square MNLS soliton fluctuations coincide completely with those in the 
case of the NLS soliton [[[], Nevertheless, there is a substantial difference between NLS and MNLS equations. 
While the correlation functions for the NLS soliton do not contain soliton parameters (in dimensionless units), 
these parameters enter explicitly into the r.h.s. of ( )l5| ) for the MNLS equation. Therefore, we can reduce the 
Gordon-Haus jitter for the MNLS soliton by varying soliton parameters. In Fig. 1 we demonstrate this effect 
in terms of variables v and w. It is seen that the function (|l^) smoothly decreases when v approaches a" 1 . 

The role of the optical energy 7 is depicted on Fig. 2. As follows from (|l7|), we have < 7 < 7r/4 for 
v < a^ 1 and 7r/4 < 7 < tt/2 for v > aT 1 . It is seen that the function ( |l6| ) grows significantly for all 7 for 
v > a -1 , while the same growth for v < a -1 is displayed to a lesser extent, including non-monotone behavior 
for small 7. 



w(l — av) = a cot 27. 



(17) 
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Figure 2: Typical behavior of the function F(a,v,j) versus v for different 7 and a = 0.05; (a) v < a , (b) 
v > a -1 . 

For slightly different definition of the soliton position as x — — — (t—q(z)) the function F(a,v,j) takes more 
simple form, 

F(a, v, 7) = 2(1 - av){l ~ cot 2 2 7 + 7 [(1 - tan 2 2 7 ) 2 + sec 2 2 7 ] cot 3 2 7 }, 

but with the same qualitative behavior. 

As regards the applicability of the above result on the soliton jitter suppression to actual femtosecond optical 
solitons, we believe that the MNLS equation is the true integrable model to start with the account for the third- 
order dispersion and the Raman self-frequency shift. It should be noted that in the adiabatic approximation 
the third-order dispersion does not contribute to the MNLS soliton velocity and width, just as for the NLS 
soliton . The Raman effect contribution can be calculated additively as an extra perturbation to the MNLS 
soliton, within the framework of the MNLS soliton perturbation theory [Q. In so doing, a value of a should 
be carefully reconciled with the specific conditions of femtosecond pulse propagation. It is worth noting that 
there exists a possibility |l7j to partially compensate for the Raman self-frequency shift effect, as applied to 
femtosecond optical pulses. Corresponding results will be published elsewhere. 

In conclusion, we have shown that for solitons of the MNLS equation which can serve as the integrable model 
for the description of ultrashort soliton dynamics, the Gordon-Haus jitter can be significantly reduced by means 
of matching soliton parameters alone. We stress that this reduction is non-perturbative with respect to a, i.e., it 
cannot be revealed in the framework of the NLS equation with a-dependent term considered as a perturbation. 
Though the above result is valid within the adiabatic approximation of the MNLS soliton perturbation theory, 
we believe it describes the main features of this phenomenon for ultrashort optical pulses. 
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